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Abstract. 

We test and compare a number of existing models predicting the location of magnetic 
reconnection at Earth’s dayside magnetopause for various solar wind conditions. We em¬ 
ploy robust image processing techniques to determine the locations where each model 
predicts reconnection to occur. The predictions are then compared to the magnetic sep¬ 
arators, the magnetic field lines separating different magnetic topologies. The predictions 
are tested in distinct liigh-resolution simulations with interplanetary magnetic field (IMF) 
clock angles ranging from 30 to 165 degrees in global magnetohydrodynamic simulations 
using the three-dimensional Block-Adaptive Tree Solarwind Roe-type Upwind Scheme 
(BATS-R-US) code with a uniform resistivity, although the described techniques can be 
generally applied to any self-consistent magnetosphere code. Additional simulations are 
carried out to test location model dependence on IMF strength and dipole tilt. We find 
that most of the models match large portions of the magnetic separators when the IMF 
has a southward component, with the models saying reconnection occurs where the lo¬ 
cal reconnection rate and reconnection outflow speed are maximized performing best. 

When the IMF has a northward component, none of the models tested faithfully map 
the entire magnetic separator, but the maximum magnetic shear model is the best at 
mapping the separator in the cusp region where reconnection has been observed. Pre¬ 
dictions for some models with northward IMF orientations improve after accounting for 
plasma flow shear parallel to the reconnecting components of the magnetic fields. Im¬ 
plications for observations are discussed. 


1. Introduction 

Magnetic reconnection occurs in plasmas where oppo¬ 
sitely directed magnetic fields merge, effectively break and 
transfer their magnetic energy into kinetic energy and heat. 
This process occurs between the interplanetary magnetic 
field (IMF) and the terrestrial magnetic field at the dayside 
magnetopause, transferring solar wind plasma into near- 
Earth space; this is the first step in the Dungey cycle of 
magnetospheric convection [ Dungey, 1961]. However, pre¬ 
dicting where magnetic reconnection occurs at Earth’s day- 
side magnetopause for arbitrary solar wind conditions has 
been the subject of studies for nearly fifty years and re¬ 
mains a challenging problem. Knowledge of where recon¬ 
nection occurs at Earth’s magnetopause is crucial for devel¬ 
oping a quantitative prediction of geomagnetic activity for 
space weather applications and to support satellite missions 
studying magnetic reconnection, such as NASA’s upcom¬ 
ing Magnetospheric Multiscale (MMS) mission [Burch and 
Drake, 2009; Moore et al, 2013]. 

In the classical sketches of Dungey [1961, 1963], recon¬ 
nection occurs at the subsolar point for due southward IMF 
and near the polar cusps for due northward IMF orienta¬ 
tions (ignoring Earth’s dipole tilt). However, much less 
is known about where reconnection occurs when the IMF 
makes an arbitrary clock angle #imf with the Earth’s dipole 
axis. The location of dayside magnetic reconnection was 
originally discussed in terms of anti-parallel [Crooker, 1979] 
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or component reconnection, either with a uniform out-of- 
plane (guide) magnetic field [Sonnerup, 1974; Gonzalez and 
Mozer, 1974] or equal and opposite components of the recon¬ 
necting magnetic field [ Cowley , 1976]); the latter is known to 
not be valid. A number of literature reviews have been writ¬ 
ten on this subject [Dorelli et al., 2007; Paschmann, 2008; 
Fuselier et al., 2011; Guo et al, 2013]; the upshot is that 
both component and anti-parallel reconnection have been 
seen in observations and simulations. This implies that nei¬ 
ther model can explain all the data so models going beyond 
them are necessary. 

Some recent models suggest that reconnection occurs 
where some parameter is maximized which produces recon¬ 
nection at an optimum efficiency. One such model is the 
maximum magnetic shear model [ Trattner et al., 2007], an 
extension of the anti-parallel reconnection hypothesis. This 
model suggests that reconnection occurs where the mag¬ 
netic shear angle between the magnetospheric and mag¬ 
netosheath magnetic fields is a maximum. Three related 
models suggest reconnection occurs where the reconnection 
outflow speed is maximized [Swisdak and Drake, 2007], the 
reconnection rate is maximized (Shay, private communica¬ 
tion, 2009; [ Borovsky, 2013]), and the reconnecting field’s 
magnetic energy is maximized [Hesse et al., 2013]. Other 
models suggest that reconnection occurs where the current 
density is a maximum [Alexeev et al., 1998] or where the 
magnetospheric magnetic field and the draped IMF are bi¬ 
sected [Moore et al., 2002]. 

There are a number of studies testing some of these mod¬ 
els in different contexts. The maximum magnetic shear 
model was shown to be consistent with observations of mag¬ 
netic reconnection at planetary magnetopauses with data 
from Cluster [Trattner et al., 2007; Dunlop et al, 2011], 
THEMIS [Trattner et al., 2012], and Cassini [Fuselier et al., 
2014]. Also, observations [Teh and Sonnerup, 2008] and sim¬ 
ulations [Schreier et al., 2010] suggest that the maximization 


1 



X - 2 


KOMAR ET AL.: TESTS OF DAYSIDE RECONNECTION MODELS 


L 



Figure 1. Schematic of an arbitrary magnetospheric 
magnetic field Bms making an angle 9 with the mag¬ 
netosheath magnetic field Bsh in the boundary normal 
(LMN) coordinate system. The angle a is the angle the 
magnetospheric field makes with the M axis. Adapted 
from Sonnerup [1974] and Swisdak and Drake [2007]. 


of the asymmetric reconnection outflow speed controls the 
location and orientation of reconnection, although it is diffi¬ 
cult to distinguish maximum outflow speed from maximum 
reconnection rate. The angle of bisection hypothesis has 
been tested with observations and in global magnetospheric 
simulations [ Pulkkinen et al., 2010] and in two-dimensional 
(2D) reconnection simulations [Hesse et al, 2013]. 

One reason the location of reconnection at the dayside 
magnetopause is difficult to predict is that reconnection does 
not occur at a single point, but rather occurs along a fine of 
finite extent. Magnetopause reconnection occurs along the 
line that separates magnetic fields of different topologies - 
solar wind and closed magnetospheric field lines merge and 
reconnect, resulting in open field lines. This line is called 
the magnetic separator [ Cowley , 1973; Stern, 1973; Siscoe, 
1987; Lau and Finn, 1990; Siscoe et al., 2001; Parnell et al., 
2010], and is the three-dimensional (3D) analog of the 2D 
X-line. There is observational [ Xiao et al., 2007; Pu et al, 
2013; Guo et al, 2013] and numerical [Hu et al, 2004; Laiti- 
nen et al., 2006, 2007; Dorelli et al., 2007; Dorelli and Bhat- 
tacharjee, 2008, 2009; Hu et al., 2009; Pulkkinen et al, 2010; 
Ouellette et al., 2010; Peng et al, 2010; Cnossen et al., 2012; 
Komar et al., 2013] evidence for the existence of separators 
throughout Earth’s magnetosphere. 

Consequently, a careful assessment of the correctness of 
the recent models for the location and orientation of recon¬ 
nection begins with a comparison of their predictions to the 
measured magnetic separators. To do so, we interpret the 
models, as Trattner et al. [2007] has done with the max¬ 
imum magnetic shear model, as predicting the collections 
of points that locally maximize a particular quantity rather 
than the single point that globally maximizes the quantity 
in question. 

The goal of this study is to systematically test the re¬ 
cent models against the separators arising in self-consistent 
global simulations of the magnetosphere. To do so, we 
employ image processing techniques to unambiguously de¬ 
termine each model’s prediction in our simulations. Each 
model’s prediction is then compared with the magnetic sep¬ 
arator determined with the algorithm described in Komar 
et al. [2013]; this algorithm has been shown to work for IMF 
with an arbitrary orientation. The present study tests each 


model’s agreement with the separators for various IMF clock 
angle #imf, IMF strength Bimf, and dipole tilt angle. We 
also present a test of one possible way to incorporate the 
effects of a flow shear parallel to the reconnecting fields into 
the models based on reconnection outflow speed and recon¬ 
nection rate. 

The layout of this paper is as follows: Section 2 reviews 
the models of dayside reconnection location that we test in 
the present study. Section 3 describes our methodology: the 
method used to determine the magnetopause, a robust im¬ 
age processing algorithm to determine each model’s predic¬ 
tion, and details of our global magnetospheric simulations. 
Our results are presented in Section 4. A summary of our 
results and applications are discussed in Section 5. 


2. Reconnection Location Models 


This section reviews the models tested in this study. 
We employ the boundary normal (LMN) coordinate sys¬ 
tem. The N direction is the magnetopause normal cor¬ 
responding to the inflow direction, the L direction cor¬ 
responds to the reconnecting component of the magnetic 
field, and the M direction gives the direction of the out-of- 
plane (guide) magnetic field, completing the right-handed 
triplet. Figure 1, adapted from Sonnerup [1974] and Swisdak 
and Drake [2007], displays the magnetospheric and magne¬ 
tosheath magnetic fields Bms and Bsh at the magnetopause 
viewed along the inflow direction. These two fields point in 
arbitrary directions in the L-M plane, making a magnetic 
shear angle of 9 with each other. The angle a is the angle 
B ms makes with the M axis and is constrained to occur 
within the range 0 < a < 9. It is tacitly assumed that the 
magnetic field components along the N axis are weak. 

In the maximum magnetic shear model [Trattner et al, 
2007], reconnection occurs at magnetopause locations where 
the magnetic shear angle 9 between the magnetospheric and 
magnetosheath magnetic fields is a maximum. This model 
identifies the anti-parallel reconnection locations of Crooker 
[1979] with high fidelity, while also identifying additional lo¬ 
cations away from the anti-parallel regions where 9 is locally 
a maximum. The magnetic shear angle 9 is calculated from 
the magnetospheric and magnetosheath magnetic fields as 


Bms • B S h 
Bms Bsh 


(1) 


Alternative explanations have suggested that reconnec¬ 
tion occurs at locations where some measure of reconnection 
efficiency is maximized. In the model of Swisdak and Drake 
[2007], reconnection on the dayside is oriented (has an angle 
a) to maximize the reconnection outflow speed CA.out■ In 
a subsequent model, reconnection orients so as to maximize 
the reconnection electric field E (Shay, private communica¬ 
tion, 2009), [ Borovsky , 2013]). To calculate the predictions 
from these models at the dayside magnetopause, one needs 
the plasma parameters of the magnetosheath and magneto¬ 
sphere, where asymmetries in the magnetic field and plasma 
density strongly affect reconnection. The scaling of asym¬ 
metric reconnection was developed for 2D antiparallel re¬ 
connection using conservation of mass and energy [Cassak 
and Shay, 2007] to obtain an expression for the asymmetric 
outflow speed 


^ A,out 


Bms,lBsh,l (Bms.l + Bsh.l) 
Po (pmsBsh.l + PshBms.l) 


and for the asymmetric electric field 


E 


CA,out 


Bms.lBsh.l 
Bms,l + Bsh.l 


26 

B’ 


( 2 ) 


(3) 
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where pms and psh are the upstream magnetospheric and 
magnetosheath plasma mass densities, respectively, and 5 
and L are the half-width and half-length of the diffusion re¬ 
gion. Equation (3) is valid independent of dissipation mech¬ 
anism [Cassak and Shay, 2008, 2009]; the present study em¬ 
ploys a resistivity p and the reconnection rate for asymmet¬ 
ric collisional, or Sweet-Parker, reconnection is [Cassak and 
Shay, 2007] 

Esp ~ ^’2“ - L?ms,l-Bsh,l. (4) 

From Fig. 1, the reconnecting component of the upstream 
magnetospheric magnetic field can be written as 

Bms,l = Sms sin a (5) 

and for the magnetosheath as 

5sh,l = Bsh sin (6 — a). (6) 

The asymmetric outflow speed is maximized by setting 

out _ 

da 



Figure 2. Magnetopause locations (green) from maxima 
in the current density’s magnitude in a simulation with 
#imf = 120°. The magnetic separator (blue) is shown for 
reference. 


and solving for a [Swisdak and Drake, 2007]. A similar ap¬ 
proach is used to maximize E (Shay, private communication, 
2009; [ Borovsky, 2013]). In practice, we interpret these as 
saying that the local CA,out and E can be found at every 
location on the magnetopause by finding a at every loca¬ 
tion, and that reconnection occurs where the parameter is a 
local maximum. We note in passing that this approach as¬ 
sumes that any nonzero Bm component of either magnetic 
held does not affect the outflow speed or reconnection rate. 
This is not likely to be correct in the real magnetosphere 
where finite Larmor radius effects are important, but should 
be accurate for the MHD simulations in the present study. 
Another limitation of the present formulation of these ar¬ 
guments is that they neglect the importance of the plasma 
bulk How parallel to the reconnecting magnetic Held com¬ 
ponents; we detail an attempt to account for this effect and 
measure the impact on the predictions of the asymmetric 
reconnection models in Section 4.4. 

More recently, Hesse et al. [2013] suggested that the re¬ 
connection rate is proportional to the product of the recon¬ 
necting component’s magnetic energies, 


E oc -Bms,lI3sh,l 5 


(7) 


vector simplifies to —r/J 2 at the reconnection site. We there¬ 
fore argue that, up to a scaling factor, the Poynting vector’s 
divergence would give similar results as the maximum cur¬ 
rent density magnitude argument for this study. 

For completeness, we also test the component reconnec¬ 
tion hypothesis as originally proposed by Sonnerup [1974] 
and Gonzalez and Mozer [1974], even though observa¬ 
tions [Gosling et al., 2007; Gosling, 2007; Teh and Sonnerup, 
2008; Phan et al., 2009, 2010; Gosling and Phan, 2013; Di- 
Braccio et al., 2013] and simulations [Schreier et al., 2010] 
suggest the model is not valid. In this model, the local recon¬ 
nection geometry orients itself such that the M component 
of the magnetic field is uniform, which, using the notation 
in Fig. 1, implies 


tan a = 


b — cos 9 
sin# 


( 8 ) 


and reconnection maximizes this product. Applying a max¬ 
imization argument with respect to a gives a solution of 
a = 9/ 2; this is equivalent to previous arguments where 
the guide field bisects the merging magnetosheath and mag¬ 
netospheric magnetic fields [Moore et al, 2002; Borovsky, 
2008]. 

We also test the model that predicts reconnection to 
occur where the current density magnitude ]J| is a maxi¬ 
mum [Alexeev et al., 1998]. A related approach has been 
to look at the divergence of the Poynting vector. Reconnec¬ 
tion has been located by measuring this quantity in observa¬ 
tions [Anekallu et al, 2013] and simulations [Papadopoulos 
et al., 1999; Palmroth et al, 2003; Laitinen et al., 2006, 2007; 
Pulkkinen et al., 2010; Palmroth et al., 2012; Hoilijoki et al., 
2014]. In steady-state, the divergence of the Poynting vector 

v C^)~- ej ' 

where E is the electric field given by Ohm’s law. The present 
study employs a large explicit resistivity, to be described in 
Section 3.4. In this system, the divergence of the Poynting 


where b = -Bms/Ssh- This model places an additional con¬ 
straint on the orientation angle a, with valid solutions only 
when the magnetic shear angle 9 is greater than the critical 
shear angle 9 C = cos -1 (1/6). 

3. Methodology 

Each of the models discussed in the previous section sug¬ 
gest reconnection at the dayside magnetopause occurs where 
some parameter is a local maximum. We model our ap¬ 
proach after Trattner et al. [2007], who identified a curve 
corresponding to where the magnetic shear angle 9 was a 
local maximum at the magnetopause. The magnetospheric 
magnetic field was modeled using the T96 model [Tsyga- 
nenko, 1995] and the Cooling model is used for the draped 
magnetosheath magnetic field [Cooling et al, 2001] to calcu¬ 
late 9 at many locations on the dayside magnetopause. The 
local shear angle 9 was plotted as a function of the mag¬ 
netopause’s y and 2 coordinates (it was projected into the 
x = 0 plane) and maxima in this shear angle map were deter¬ 
mined by finding maxima along cuts parallel to the IMF (K. 
Trattner, private communication, 2013). 

The following sections detail how we robustly identify the 
magnetopause and calculate model predictions in our global 
magnetospheric simulations. 

3.1. Identification of the Magnetopause and Plasma 
Parameter Sampling 
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To identify the magnetopause in our global magneto- 
spheric simulations, we modify the method of Nemecek 
et al. [2011]. (One may alternatively adopt the method 
described in Palmroth et al. [2003] to identify the magne¬ 
topause using streamlines, although this is not tested here.) 
The current density magnitude is sampled radially from 
7 < r < 20 Earth radii (Re) at our highest grid cell resolu¬ 
tion (0.125 Re). The location of maximum current density 
magnitude at smallest r is identified as the magnetopause 
and this location’s coordinates tmp are saved; the chosen 
sampling range excludes the ring current, but we take ad¬ 
ditional care by ensuring each point has a radial distance 
within 2 Re of the previous point’s radial distance. We 
employ spherical coordinates in the Geocentric Solar Mag- 
netospheric (GSM) system, where the azimuthal angle <j> is 
measured from the +x axis and the polar angle 9 is mea¬ 
sured from the +z axis; we discretize both into 5° incre¬ 
ments. (Other increments were tested and found to con¬ 
verge to the 5° results.) We use an azimuthal range of 
— 130° < <j> < 130° to map the magnetopause for northward 
IMF orientations (Simp = [30°, 60°, 90°] for the purposes of 
this study) and —110° < (j> < 110° for southward IMF ori¬ 
entations (0imf = [120°, 150°, 165°]); the polar angle has 
a range 0° < 9 < 180° for both northward and southward 
IMF orientations. 

An example of this process can be seen in Fig. 2, display¬ 
ing the magnetopause surface as green spheres for a sim¬ 
ulation with IMF clock angle 120° that will be described 
in Section 3.4. Additionally, the separator for this simu¬ 
lation from Komar et al. [2013] is plotted in blue. It lies 
within the magnetopause surface, as it should, confirming 
this method. The grid is rotated out of the ecliptic plane; 
the spherical coordinate system is rotated counter-clockwise 
about the +x axis by an angle 0imf/2. This is done to 
achieve near uniform sampling near the separator and each 
model’s prediction, i.e. near 9 « 90° in the rotated frame. 

Once the approximate location of the magnetopause has 
been determined, the local plasma parameters of the magne¬ 
tosphere and magnetosheath are sampled. The appropriate 
way to do this is to sample along the magnetopause nor¬ 
mal n at each magnetopause position tmp. One may per¬ 
form a minimum variance analysis (MVA) on the magnetic 
field [Sonnerup and Cahill, 1967] to determine the magne¬ 
topause normal. However, we find in practice this results in 
incorrect magnetopause normals since the magnetospheric 
and magnetosheath magnetic fields are parallel near the sep¬ 
arator, violating the core assumption of MVA. 

Instead, the normal is calculated with the method de¬ 
scribed in Hoppe et al. [1992]. The algorithm is to take a 
single point on the magnetopause tmp and find its four near¬ 
est neighbors r*, on the magnetopause, with k £ [1, 2, 3, 4]. 
Calculate the difference vector r*, — tmp, which gives the 
distance from the point on the magnetopause and its fcth 
nearest neighbor. Then, construct the covariant matrix M, 
whose ij element is 

Mij (r M p) = T, (r fe - r M p) ; (r fe - rMp)^ , (9) 

k 

where i and j refer to the x, y, or 2 component in GSM co¬ 
ordinates. The eigenvector corresponding to the minimum 
eigenvalue of M is the magnetopause normal; this normal 
direction is chosen to point away from Earth, i.e., it - n > 0, 
and has been confirmed by inspection. This procedure has 
a physical interpretation of placing several planes at tmp 
with arbitrary orientations; the normal of the plane that 
minimizes the distance between the plane and the nearest 
neighbors is the magnetopause’s normal. 

The plasma parameters are sampled along the normal 
vector n at our highest resolution (0.125 Re) between 
tmp — (5 Re) n and tmp + (5 Re) n. To accurately measure 



Figure 3. Magnetic shear angle 9 (blue spheres) as a 
function of the magnetopause’s y and 2 coordinates and 
the corresponding ridge of maximum magnetic shear an¬ 
gle (red spheres) in a simulation with #imf = 120°. 

the asymptotic plasma densities and magnetic fields under¬ 
going magnetic reconnection, we employ the method used 
by Cassak and Shay [2009] for systems with asymmetries in 
either parameter. The current density is measured along n 
to determine where it falls to 1/e of its maximum value in 
the magnetosheath. The distance between the current maxi¬ 
mum and this location is the current sheet’s half-thickness 5. 
The upstream plasma parameters of the magnetosphere are 
measured at tmp — (2<5) n and those of the magnetosheath 
at tmp + (2 5) n. We confirm the measured magnetospheric 
and magnetosheath magnetic fields and densities are rep¬ 
resentative of the asymptotic values, as desired. The raw 
upstream magnetic fields can have a small normal compo¬ 
nent (\B n \ ~ 1 nT); these are removed prior to calculations. 
These values are used as Bms, Bsh, Pms, and psh for the 
models discussed in Section 2. 

We note in passing that identifying the magnetopause 
as the maximum in the current density’s magnitude has 
its limitations. One consequence of asymmetric reconnec¬ 
tion systems is that the stagnation point separates from the 
reconnection X-line in 2D reconnection [Cassak and Shay, 
2007]; the maximum of the out-of-plane current density also 
separates from the X-line. We will describe in Section 4.2 a 
test of reconnection location models in distinct simulations 
with IMF strengths of 20, 5, and 2 nT. The separation dis¬ 
tances between the current density maximum and separator 
(the 3D version of the X-line) at the subsolar point in these 
simulations increase as the IMF strength is lowered, with 
distances of 0.07 (20 nT), 0.36 (5 nT), and 0.51 Re (2 nT); 
the separation distances are resolved in the 5 and 2 nT sim¬ 
ulations since our maximum resolution is 0.125 Re- In such 
systems, the current density maximum does not coincide 
with the magnetopause. For the purposes of this study, 
we assume this feature corresponds to a compression of the 
magnetopause in the n direction at every point with minimal 
impact to each model’s calculation as the sampling method 
described here measures the plasma parameters far enough 
away from the magnetopause location. 

3.2. Finding Model Predictions Using Image Processing 
Techniques 

Having found the approximate location of the magne¬ 
topause and the parameters locally governing reconnection 
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(pMS, Psh, Bms, Bsh), the relevant quantities for each 
model are calculated at every point on the magnetopause. 
The magnetic shear angle 9 is calculated from Eq. (1) at 
every point and is retained. Then at every point, a is dis¬ 
cretized in the range 0 < a < 9. For each value of a, 
Bms,l and Bsh,l are calculated from Eqs. (5) and (6) , 
and CA,out, Esp, and Bms,lBs H , l are calculated for that 
a. The a which maximizes the quantity in question is re¬ 
tained, and the value of CA,out, Esp, and B^q h Bg H L is 
retained for that point. (One could also determine a from 
9ca, out/da = 0 and dE/da = 0.) The current density mag¬ 
nitude is also retained at every point. 

The results for each model are projected into the x = 0 
plane, resulting in a 2D image y (y, z), where J? represents 
the quantity in question for each model. An example is 
shown in Fig. 3 for the maximum magnetic shear model, 
where J? corresponds to 9, which is plotted as the third 
dimension in blue spheres for the #imf = 120° simulation. 
The magnetic shear angle surface has a saddle shape, as 
do the surfaces for all the other models tested (not shown). 
This plot makes it clear that finding the collection of loca¬ 
tions where the quantity in question is locally a maximum 
is tantamount to finding the ridge of local maxima in y. 

We achieve this in a robust way using image process¬ 
ing techniques using a ridge finding algorithm [Linde- 
berg, 1993, 1998]. First, we construct the Hessian tensor 
Hij ( y, z) = d 2 y (y, z) /dxidxj via finite differencing at ev¬ 
ery point in the image. Diagonalizing the Hessian matrix 
defines a ( p, q) coordinate system, where the eigenvector as¬ 
sociated with the maximum absolute eigenvalue points in 
the direction of the local maxima. (Without loss of gen¬ 
erality, we designate that eigenvector as p, the unit vector 
for the p direction.) Local maxima are determined by find¬ 
ing all locations where dy/dp = 0 (local extrema test) and 
d 2 y/dp 2 < 0 (local maxima test) using linear interpolation 
between nearest neighbors [ Lindeberg , 1993]. 

The results of this technique for the maximum magnetic 
shear model for 0imf = 120° are shown as the red spheres 
in Fig. 3. These red spheres clearly mark the ridge of the 
magnetic shear surface as desired. The methodology we em¬ 
ploy yields qualitatively similar results to previous studies 
testing the maximum magnetic shear model for compara¬ 
ble magnetospheric and solar wind conditions (see Fig. 1 
in Fuselier et al. [2011]), indicating this is an appropriate 
technique to determine each model’s prediction. A more 
thorough comparison between these techniques will be pre¬ 
sented in Section 4.1. 

3.3. Determining Magnetic Separators 

We employ the separator mapping algorithm of Komar 
et al. [2013] which has been shown to reliably trace the day- 
side portion of the magnetic separators connecting the mag¬ 
netic nulls in global magnetosphere simulations for any IMF 
direction. When the IMF has a northward component, ob¬ 
servations have shown that reconnection occurs poleward 
of the magnetic cusps (see Wilder et al. [2014], and ref¬ 
erences therein). We therefore trace portions of the mag¬ 
netic separator nightward of the nulls in our simulations 
with #imf = 30°, 60°, and 90°. 

In the separator tracing algorithm, a hemisphere is ini¬ 
tially centered around a magnetic null. The hemisphere’s 
surface, of radius 1 Re for our purposes, is discretized into 
a grid. The magnetic field lines piercing the hemisphere 
at each grid point are traced to determine their magnetic 
topology: closed terrestrial, solar wind, and open either 
connected to the north or south magnetic poles. The ap¬ 
proximate location of the separator is determined by finding 
where these four magnetic topologies meet on a hemisphere’s 
surface. Then, another hemisphere is centered at the deter¬ 
mined separator location, and the procedure is iterated to 
trace the separator. The dayside separator is traced from 



Figure 4. Cross polar cap potential difference A$pc 
in simulations with |Bimf| = 2, 5 and 20 nT for 
Simf = 120°. The potential difference is linear in ]Bimf|, 
so the polar cap is not saturated. 


northern to southern null in this fashion and the nightside 
portion is traced back for 10 Re from the magnetic nulls 
when the IMF has a northward orientation. 

3.4. Magnetospheric Simulation Study 

We perform global simulations using the Block Adap¬ 
tive Tree Solarwind Roe-type Upwind Scheme (BATS-R-US) 
code [Powell et al., 1999; Gombosi et al., 2000; De Zeeuw 
et al., 2000; Toth et al., 2012], although the methods detailed 
in the previous sections can be adapted to other global mag¬ 
netospheric codes. Simulations are performed at NASA’s 
Community Coordinated Modeling Center (CCMC) and are 
partially analyzed with CCMC’s Kameleon software suite. 
The simulation domain is —255 < x < 33, —48 < y < 48, 
and —48 < z < 48, where distances are measured in Re and 
the coordinate system is GSM. 

The simulations are carried out using BATS-R-US version 
8.01. The simulations are evolved for two hours (02:00:00) 
of magnetospheric time. We look at the 02:00:00 mark of 
simulation data as the system has achieved a quasi-steady 
state by this time. The standard high-resolution grid for 
CCMC simulations has 1,958,688 grid cells with a coarse 
resolution of 8 Re in the far magnetotail, and a fine resolu¬ 
tion of 0.25 Re near the magnetopause. The present study 
employs a higher resolution grid of 0.125 Re packed in the 
region —15 < x, y, z < 15 Re, totaling 16,286,400 simula¬ 
tion grid cells. 

The initial simulations do not employ a dipole tilt, al¬ 
though we include one later. All simulations use fixed so¬ 
lar wind inflow conditions. The solar wind has tempera¬ 
ture T = 232,100 K, IMF strength 20 nT, number density 
n = 20 cm -3 , and a solar wind velocity of v = —400 km/s x. 
We perform distinct simulations with IMF clock angles 
0 imf = 30°, 60°, 90°, 120°, 150°, and 165°. The IMF does 
not have a B x component. Constant Pederson and Hall con¬ 
ductances of 5 mhos are used. The solar radio flux F10.7 
index is set at a value of 150. 

The parameter regime of the solar wind in the present 
global magnetospheric simulations employed atypical IMF 
strengths and solar wind densities; both quantities are an 
order of magnitude higher than typically quiet solar wind 
conditions. This choice is made so the dayside magneto¬ 
sphere is smaller, and the region of high resolution need not 
be as large. To ensure that the results are reliable, it is im¬ 
portant to check that the polar cap is not saturated with 
these solar wind parameters [Lopez et al, 2010]. 

We measure the polar cap potential difference in three 
distinct simulations with Bimf = 2, 5, and 20 nT. These 
simulations have no dipole tilt, the IMF clock angle is fixed 
at 120°, and all other solar wind parameters are the same 
as before. Figure 4 displays the measured polar cap po¬ 
tential difference as a function of ]Bimf|. The plot clearly 
reveals that the potential difference remains linear as a func¬ 
tion of IMF strength, implying the polar cap is not satu¬ 
rated. For comparison, Lopez et al. [2010], using LFM sim¬ 
ulations, showed that the polar cap potential difference re¬ 
mains linear for Bimf < 11.25 nT with solar wind density 
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Figure 5. Results for the maximum magnetic shear model in simulations with #imf of (a) 30°, (b) 60°, 
(c) 90°, (d) 120°, (e) 150°, and (f) 165°. The calculated magnetic shear angle 6 is the color background, 
the gray squares display the curve of maximum magnetic shear angle, and the magnetic separator is 
displayed in white. The gray oval displays the magnetopause’s projection in the x = 0 plane. Maxima 
determined from cuts along the IMF direction [Trattner et al, 2007] are displayed as black squares in 
panel (f). 
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Figure 6. Results for the asymmetric outflow speed CA,out ■ See the caption of Figure 5 for definitions. 


nsw = 5 cm ”' 5 and for Bimf < 15 nT when nsw = 8 cm ” 3 
using a solar wind speed of 400 km/s and ionospheric Ped¬ 
ersen conductance of 5 mhos. As our simulations employ 
nsw = 20 cm” 3 , it is reasonable that we find the polar cap 
is not saturated for 5 imf = 20 nT. Our measured viscous 
potential (the y-intercept) is 43.8 kV, consistent with the 
trend seen by Lopez et al. [2010], with their measured vis¬ 


cous potentials of 33.9 and 41.9 kV for nsw = 5 and 8 cm” 3 , 
respectively. We are therefore confident that the polar cap 
is not saturated for our simulation parameters. 

For the present simulations, we employ a uniform ex¬ 
plicit resistivity 77 . It is known that Earth’s magnetopause 
is collisionless, but including an explicit resistivity allows for 
reproducible results that are independent of the numerics. 
With the exception of the asymmetric Sweet-Parker recon- 
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Figure 7. Results for the asymmetric Sweet-Parker reconnection rate Esp ■ See the caption of Figure 5 for definitions. 



Figure 8 . Results for the angle of bisection. See the caption of Figure 5 for definitions. 


nection rate [Eq. (4)], all of the models selected for test¬ 
ing are independent of the reconnection dissipation mech¬ 
anism. Therefore, the model predictions should likewise 
remain independent of this choice for any self-consistent 
magnetospheric model. We include an explicit resistivity 
rj/fj .0 = 6.0 x lO 10 m 2 /s in our simulations, which allows us 
to sufficiently resolve the dayside magnetopause to deter¬ 
mine the validity of the models tested. We refer the reader 
to Komar et al. [2013] for a more detailed discussion on the 
inclusion of an explicit resistivity. 


3.5. Simulation Resolution 

The numerical grid described in the present study is dif¬ 
ferent from those of Komar et al. [2013] where the algo¬ 
rithm for finding separators was presented. The present 
simulations employed a resolution of 0.125 Re in the re¬ 
gion —15 < x, y, z < 15 Re whereas Komar et al. [2013] 
used a resolution of 0.125 Re in the region —6 < x < 10, 
—10 < y < 10, —5 < z < 5 Re, totaling a difference of 
~12.6 million simulation grid cells. The present study uses 
the larger high-resolution domain to ensure that the entire 
dayside magnetopause (and therefore the magnetic separa¬ 
tors) is within the high-resolution domain. The magnetic 
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separators in the present simulations are not significantly 
different than those in Komar et al. [2013]. 

4. Results 

We present results for the predicted locations for the mod¬ 
els discussed in Section 2 and compare the results to mag¬ 
netic separators. Figures 5-9 display the results in a common 
format. The quantity of interest is displayed as the back¬ 
ground color, the model’s prediction is displayed with solid 
gray squares, and the magnetic separators determined with 
the algorithm described in Komar et al. [2013] are displayed 
as the solid white line. The magnetopause’s intersection 
with the x = 0 plane is displayed as the gray oval. The 
panels for each figure correspond to IMF clock angles #imf 
of (a) 30°, (b) 60°, (c) 90°, (d) 120°, (e) 150°, and (f) 165°. 
Figure 5 displays the magnetic shear angle 9 in degrees cal¬ 
culated from Eq. (1), Fig. 6 has the asymmetric outflow 
speed CA.out in km/s calculated from Eq. (2), Fig. 7 has the 
asymmetric Sweet-Parker reconnection rate Esp in mV/m 
calculated from Eq. (4), Fig. 8 has the quantity -Bms,l£>!h,l 
from Eq. (7) in 10 s nT 4 , and Figure 9 has the current den¬ 
sity magnitude |J| in nA/m 2 . 

Figure 10 displays the results for the uniform Bm com¬ 
ponent model. Plots (a) through (f) are for the same IMF 
clock angle simulations as previously discussed, with a few 
notable differences. The local Bm component (guide field) 
is calculated by determining the orientation angle a from 
Eq. (8), and rotating about the determined boundary nor¬ 
mal n. The gray arrows in Fig. 10 display the projection of 
the magnetic field’s M component in the x — 0 plane; loca¬ 
tions without any vector displayed are where reconnection 
is predicted to be geometrically impossible. The magnetic 
separator is displayed as the solid black line, with the mag¬ 
netopause’s location in the x = 0 plane plotted in gray as 
before. 

Section 4.1 discusses how the predictions of each model 
change with IMF clock angle for fixed solar wind conditions. 
We test the dependence on IMF strength of various models 
in Section 4.2. Section 4.3 describes results from two dis¬ 
tinct simulations that include a dipole tilt with IMF clock 
angles 120° and 30°. Lastly, Section 4.4 discusses an at¬ 
tempt to incorporate the effects of a plasma flow shear on 
reconnection into the outflow speed and reconnection rate 
models from Section 2. 

4.1. Clock Angle Dependence 

Before assessing each model, we perform a more care¬ 
ful comparison between our ridge finding technique and the 
one in Trattner et al. [2007]. We plot the result using the 
technique of Trattner et al. [2007], which finds maxima of 
magnetic shear angle along cuts parallel to the IMF direc¬ 
tion, as black squares in Fig. 5(f) for the #imf = 165° 
simulation. While qualitatively similar, there are appar¬ 
ent differences. If one takes cuts along different directions 
(not shown), the detected locations change, and can move 
to locations on the other side of the image processing ridge. 
This underscores two points: (1) it is important to employ 
a user-independent technique, such as image processing, to 
locate each models’ prediction as the Trattner et al. [2007] 
algorithm exhibits sensitivity to the cut orientation and (2) 
the noticeable changes in the prediction are caused by the 
relative flatness in the saddle region that occurs near the 
subsolar point in these results. Consequently, there are un¬ 
certainties that could cause minor changes to the predictions 
which should be kept in mind in the interpretation of these 
results. 

We now focus on the predictions of each model and how 
they compare with the separators for all IMF clock an¬ 
gles. Of the models tested, the maximum magnetic shear 
as shown has the largest discrepancy with the dayside por¬ 
tion of the magnetic separators in simulations without a 


dipole tilt, but could have better agreement in light of 
the uncertainties discussed above. More importantly, only 
the maximum magnetic shear model reproduces the portion 
of the magnetic separator nightward of the magnetic nulls 
with high fidelity. These portions of the magnetic separa¬ 
tors coincide with the antiparallel regions and are near the 
magnetic cusps, consistent with early expectations [ Dungey , 
1961, 1963; Crooker, 1979] and with observations (see Wilder 
et al. [2014], and references therein). 

The asymmetric reconnection arguments ( CA,out and 
Esp) show better agreement with the magnetic separators 
overall for the simulations with southward IMF orientations. 
The ridges of these models rotate around the x axis, as 
does the separator, although there are differences between 
the separators and the model predictions for clock angles 
#imf V 120°. 

The predictions for the angle of bisection and the max¬ 
imum current density models have a similar shape as the 
other models’ predictions through the subsolar point. Fig¬ 
ures 8(a)-(c) and 9(a)-(c) show intersections between the 
model predictions and magnetic separators near the antipar¬ 
allel regions and at the subsolar point for 9 imf < 90°. The 
predictions of these arguments do follow the separator more 
closely for southward IMF orientations, and do not rotate 
around the magnetopause for northward orientations. 

The uniform Bm component model would accurately pre¬ 
dict the separator if the arrows point along the separator. 
They clearly do not; at the subsolar point, the measured 
orientation angle a relative to the magnetospheric magnetic 
field is (see Fig. 1): (a) N/A (reconnection is geometrically 
impossible), (b) 46.9°, (c) 53.1°, (d) 64.1°, (e) 76.8°, and (f) 
83.6°. These angles are inconsistent with the orientation an¬ 
gles asep of the magnetic separators (a) 10.3°, (b) 21.7°, (c) 
35.8°, (d) 55.3°, (e) 71.2°, and (f) 77.5°. The discrepancies 
between the magnetic field’s M component and the orienta¬ 
tion of the magnetic separators provide further evidence that 
this description does not accurately predict the orientation 
of magnetic reconnection, consistent with previous observa¬ 
tions [Teh and Sonnerup, 2008] and simulations [Schreier 
et al., 2010]. 

4.2. IMF Strength Dependence 

Equations (2) and (4) have strong dependencies on 
the magnetospheric and magnetosheath magnetic field 
strengths. The magnetospheric magnetic field strength is 
usually much larger than that of the magnetosheath, and 
does not vary much for the chosen simulation parameters. 
This implies that the reconnection efficiency models depend 
strongly on the magnetosheath magnetic field strength Bsh- 
Meanwhile, the maximum magnetic shear model is expected 
to be independent of Bsh- We therefore vary £>sh to see if 
it helps to distinguish among these models. 

We lower the IMF strength from 20 nT to 5 and 2 nT 
in separate simulations and use #imf = 120° with all other 
simulation parameters unchanged. Figure 11 displays the 
results from these simulations. The top row displays the 
magnetic shear angle in degrees, the middle row has the 
asymmetric outflow speed in km/s, and the bottom row has 
the asymmetric Sweet-Parker reconnection rate in mV/m. 
Figures 11(a), (d), and (g) are for the Bimf = 20 nT sim¬ 
ulation, (b), (e), and (h) have Bimf = 5 nT, and (c), (f), 
and (i) use Bimf = 2 nT. As before, each model’s prediction 
is plotted as the gray squares, the magnetic separators are 
displayed in white. 

We note in passing that the magnetopause is elongated 
along the IMF direction in the 20 nT plots of Fig. 11. This 
effect is consistent with observations [Lavraud et al., 2013] 
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Figure 9. Results for the current density magnitude |J|. See the caption of Figure 5 for dehnitions. 
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Figure 10. Results for the uniform Bm component model. Panels (a) through (f) are from the same 
simulations as in Figures 5-9. The y and z components of Bm are plotted as gray arrows and the magnetic 
separator is displayed in black. 


and simulations [Lavraud and Borovsky, 2008] for low solar 
wind Alfven Mach number. 

The subsolar region in each simulation becomes increas¬ 
ingly asymmetric as the IMF strength is lowered; the mea¬ 
sured magnetic field ratios are b = Bms/Bsh = 1.4, 4.1, 
and 11 for Bimf = 20, 5, and 2 nT, respectively. As seen in 
Fig. 11, each model’s prediction becomes noisier for smaller 
Bsh due to the limitation of determining the magnetopause 
via the current maximum, as described in section 3.1. 


Figures ll(a)-(c) show the maximum magnetic shear 
model remains relatively independent of the IMF strength as 
expected; this model’s predictions have very similar struc¬ 
tures as seen previously for arbitrary #imf, making simi¬ 
lar angles with the magnetospheric magnetic field of 49.6°, 
50.2°, and 53.2° in simulations with Bimf = 20, 5, and 
2 nT, respectively. This model maps large portions of the 
magnetic separator across the dayside magnetopause, con¬ 
sistent with previous results. 

It is difficult to see differences between the asymmetric 
outflow speed and asymmetric Sweet-Parker reconnection 
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Figure 11. Results for testing model dependence on IMF strength for: (top row) the magnetic shear 
angle, (middle row) asymmetric outflow speed, and (bottom row) the asymmetric Sweet-Parker recon¬ 
nection rate. Plots (a), (d) and (g) are for Bimf = 20 nT, (b), (e), (h) are for Bimf = 5 nT, and (c), (f), 
(i) for Bimf = 2 nT. See the caption of Figure 5 for definitions. 


rate models in Fig. 11. The predictions of both models 
closely map the magnetic separator in these simulations. 
For a more quantitative comparison between these models, 
we calculate each model’s predicted orientation angle and 
compare these values to the separator’s orientation at the 
subsolar point. The separator’s orientation angle as ep mea¬ 
sured from the magnetospheric magnetic field is 55.3°, 61.5°, 
and 56.2° in the Bimf = 20, 5, and 2 nT simulations, respec¬ 
tively. The orientation of the ridges passing through the sub¬ 
solar region for the asymmetric outflow speed make an angle 
ay and, likewise, the asymmetric Sweet-Parker reconnection 
rate ole are again measured from the magnetospheric mag¬ 
netic field, with values ay = 47.0° and Oe = 51.0° in the 
Bimf = 20 nT simulation, ay = 45.7° and a® = 49.7° for 
Bimf = 5 nT, and finally, ay = 47.4° and as = 47.4° with 
Bimf = 2 nT. 

The orientation angle as decreases as IMF strength is 
lowered, while ay displays a minimum when Rimf = 5 nT. 
However, asep is maximum when £?imf = 5 nT. All of the 
models tested have predictions that are reasonably close, 
but none show perfect agreement with the orientation at 
the subsolar point. This implies that these models will need 
further modification to accurately predict the orientation of 
reconnection at Earth’s subsolar magnetopause. 

4.3. Dipole Tilt Dependence 


Including a dipole tilt can alter the location of magnetic 
reconnection. For example, Trattner et al. [2007] showed 
that the prediction of the maximum magnetic shear model 
moves southward under positive dipole tilt (as measured 
from the +z axis in the x-z GSM plane, i.e. the north¬ 
ern magnetic pole pointing Sunward), and likewise moves 
northward for negative dipole tilt when paired with south¬ 
ward IMF orientations. Reconnection in simulations has 
also been measured to change location in simulations with a 
dipole tilt [Park et al., 2006; Palmroth et al., 2012; Cnossen 
et al, 2012; Hoilijoki et al, 2014]. 

Employing a more realistic dipole model of Earth may as¬ 
sist in distinguishing among the models. Note in Figs. 5-9 
the model predictions intersect the magnetic separators at 
the subsolar point and near the anti-parallel region for all 
IMF clock angles. These intersections result from symme¬ 
try; including a dipole tilt will break this symmetry making 
it easier to distinguish which model predicts the magnetic 
separator. 

We present the results of two additional simulations that 
employ a +15° dipole tilt. The same solar wind parameters 
described in Section 3.4 are used in separate simulations 
with IMF clock angles of 120° (Fig. 12) and 30° (Fig. 13). 
Figures 12 and 13 display each model’s prediction as the grey 
squares, where the maximized quantity is: (a) the magnetic 
shear angle 9 in degrees, (b) the asymmetric outflow speed 
CA,out in km/s, (c) the asymmetric Sweet-Parker reconnec¬ 
tion rate Esp in mV/m, (d) the quantity B^ IS L 1T? H l from 
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Figure 12. Results from a simulation with a dipole tilt of +15° and a southward IMF orientation with 
clock angle 120°: (a) the magnetic shear angle 9, (b) the asymmetric outflow speed CA,out, (c) the asym¬ 
metric Sweet-Parker reconnection rate Esp, (d) the angle of bisection, (e) the current density magnitude 
| J|, and (e) the y and 2 components of the uniform Bm component. See the caption of Figure 5 for 
definitions. 
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Figure 13. Results from a simulation with dipole tilt of +15° and a northward IMF orientation with 
clock angle 30°. See the caption of Figure 12 for definitions. 


Eq. (7) in 10 8 nT 4 , and (e) the current density magnitude 
|J| in nA/m 2 , and (f) the projection of the uniform Bm 
component in the x = 0 plane. The magnetic separators are 
displayed in white for (a)-(e) and black in (f) as before. 

In Fig. 12, each model appears to map the magnetic sepa¬ 
rator reasonably well with a southward IMF orientation and 


positive dipole tilt. The maximum magnetic shear model’s 
prediction follows the separator more closely than any other 
model when the dipole tilt is present. The other models map 
the separator more closely in the northern, dusk quadrant 
in Figs. 12(b)-(e) than any other quadrant. [We use dawn 
(— y), dusk (+t/), southern (— 2 ),, northern (+ 2 ) to indicate 
position in the y-z GSM plane.] Figure 12(f) shows that 
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the uniform Bm component model again fails to map the 
magnetic separator in any quadrant. 

The results for IMF clock angle 30° and positive dipole 
tilt in Fig. 13 are significantly different. The separa¬ 
tor moves approximately 45° in longitude duskward of the 
+2?gsm axis. The model predictions intersect the mag¬ 
netic separator in the northern, dusk and southern, dawn 
quadrants, with the maximum current density mapping the 
largest portion of the separator in the northern, dusk quad¬ 
rant. The maximum magnetic shear model is again the only 
model that faithfully reproduces the nightward portion of 
the magnetic separator. Most importantly, no model maps 
any part of the separator in the southern, dusk quadrant. 

4.4. Effect of a Flow Shear 

Existing models have thus far ignored the impact of solar 
wind flow around the magnetopause. The tailward flow of 
magnetosheath plasma around the magnetopause can have 
a component parallel to the reconnecting component of the 
magnetic field, which is known to impact the efficiency of 
reconnection [Chen and Morrison, 1990; La Belle-Hamer 
et al, 1994, 1995; Cassak and Otto, 2011]. This effect po¬ 
tentially could alter models based on reconnection efficiency. 

The incorporation of a plasma flow shear into equa¬ 
tions (2) and (4) is ongoing work. Based on preliminary 
work (C. E. Doss et al., private communication, 2014), we 
hypothesize that the predictions for the outflow speed [Cas¬ 
sak, 2011] and reconnection rate [Cassak and Otto, 2011] 
for symmetric reconnection with a flow shear can be gen¬ 
eralized for asymmetric reconnection by replacing the sym¬ 
metric outflow speed ca with the asymmetric outflow speed 
CA,out from Eq. (2), giving 


CShearflow — CA,out 




n“ 

° A,out 


( 10 ) 


shear, and remain largely unchanged with a southward IMF 
orientation. Quantifying this effect in the northward IMF 
simulation, the prediction for CA,out makes an angle a of 
44.7° with the magnetospheric magnetic field, and changes 
to 32.7° after accounting for the flow shear; the separator 
makes an angle of as ep = 10.3°. Likewise, the prediction 
of Esp makes an angle of 43.2° and is 29.7° after includ¬ 
ing a flow shear. Changes to the orientation angles for 
each model’s prediction in the southward IMF simulation 
are much smaller changing from 47.0° to 46.4° for the out¬ 
flow speed, and change from 51.0° to 52.1° for the reconnec¬ 
tion electric field; the separator has asep = 55.3° for this 
simulation. A further test of this model was performed for 
the 30° clock angle with positive dipole tilt simulation. The 
results (not shown) revealed improvement of the model pre¬ 
dictions in some areas, and worse agreement in others, but 
still differed significantly from the separator. 

The relative importance of a flow shear is understood as 
a result of general flow patterns around the magnetopause. 
For southward IMF, reconnection occurs near the ecliptic 
plane where the magnetosheath flow and magnetospheric 
convection are tangential to this plane. Thus, any flow 
shear parallel to the magnetic field is small and reconnec¬ 
tion should not be strongly affected. However, the magne¬ 
tosheath flow at high latitudes is quite large. When the IMF 
is northward, there is appreciable flow parallel to the recon¬ 
necting magnetic fields and may lead to the suppression of 
reconnection. The present results suggest that including 
flow shear improves the agreement between separators and 
model predictions for the reconnection efficiency arguments. 

We note that the analysis here assumes that any vm flow 
tangential to the L-components of the magnetic fields has 
no effect on the reconnection. This is unlikely to be the case, 
but including this is beyond the scope of the present study. 

5. Conclusions 


and 

Eshearflow — EsP [ 1 ~n 

V C A,out 

where Esp is given by Eq. (4). The flow shear speed v s is 
given by 

vsh.l — ums.l 

v a = - 1 -— 

2 

and ush.l and ums.l are the magnetosheath and magne¬ 
topause flow in the L direction, respectively. 

We test to see if these expression alter the predictions of 
the reconnection efficiency arguments. As before, a is dis¬ 
cretized to calculate the asymmetric arguments CA.out and 
Esp while the L direction is rotated around the magne¬ 
topause normal for each a. The velocities on each side of 
the magnetopause are measured 3<5 upstream along n and 
vsh,l and ums.l are calculated by taking the L-component 
for each a. We calculate the modified reconnection outflow 
speed and electric field using equations (10) and (11) as a 
function of a. (For v s > c A ,out, eshearflow and EshearfW are 
set to zero as reconnection should not be possible.) Modi¬ 
fied model predictions are determined with image processing 
techniques as described in Section 3.2. 

Figure 14 displays results for two different clock angles: 
the top row are from the simulation with #imf = 30° and 
the bottom row has 0 imf = 120°. Panels (a) and (e) are 
results without corrections for flow shear and are replicated 
from Fig. 6. Likewise, panels (c) and (g) are from Fig. 7 for 
the reconnection rate. Panels (b) and (f) give results that 
incorporate a flow shear in the calculation of the outflow 
speed and (d) and (h) on the reconnection rate. The model 
predictions more closely map the separator when the IMF 
has a northward orientation after accounting for the flow 



In this paper, we calculate the predictions of several mag¬ 
netic reconnection location models at the dayside magne¬ 
topause in global resistive MHD simulations. Models tested 
are: the maximum magnetic shear model [Trattner et al., 
2007], maximization of the asymmetric reconnection out¬ 
flow speed [Swisdak and Drake, 2007], maximization of the 
asymmetric reconnection rate (Shay, private communica¬ 
tion, 2009; [Borovsky, 2013]), the angle of bisection [Moore 
et al., 2002; Borovsky, 2008; Hesse et al, 2013], the max¬ 
imization of the current density magnitude [Alexeev et al., 
1998], and component reconnection (uniform Bm compo¬ 
nent) [Sonnerup, 1974; Gonzalez and Mozer, 1974]. The 
first five models predict reconnection to occur at locations 
where a given parameter is maximized. We employ robust 
image processing techniques to determine each model’s pre¬ 
diction. Each model is tested in separate global magne¬ 
tospheric simulations with various IMF orientations, IMF 
strengths, and with a dipole tilt. The predictions for each 
model are then compared to the magnetic separators to de¬ 
termine which model, if any, accurately predicts the location 
of reconnection at the dayside magnetopause. We summa¬ 
rize our findings: 

1. For the chosen simulation parameters, we find that all 
models are within a few Earth radii (Re) of the magnetic 
separators when the IMF has a southward orientation and 
no dipole tilt. 

2. None of the models chosen for this study faithfully re¬ 
produce the entire magnetic separators when the IMF has 
a northward orientation and no dipole tilt. However, only 
the maximum magnetic shear model faithfully reproduces 
the portion of the separator nightward of the magnetic nulls 
where reconnection is expected to occur. 
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Figure 14. Results for testing the effects of flow shear on model predictions with (top row) #imf = 30° 
and (bottom row) 120°. Plots (a) and (e) are for the unmodified asymmetric outflow speed CA,out, (b) 
and (f) are for the modified outflow speed cshearflow, (c) and (g) are for the unmodified asymmetric 
Sweet-Parker reconnection rate Esp, and (d) and (h) are for the modified reconnection rate Eshearflow- 
See the caption of Figure 5 for definitions. 


3. The asymmetric outflow speed and asymmetric Sweet- 
Parker reconnection rate more closely map the magnetic sep¬ 
arators than other models tested under southward IMF ori¬ 
entations and no dipole tilt. The maximum magnetic shear 
model has a fixed orientation at the subsolar magnetopause 
and does not rotate with the magnetic separators for differ¬ 
ent IMF clock angles in simulations without a dipole tilt. 

4. As the IMF strength is varied with fixed clock angle, 
the magnetic separators display different orientations than 
the predictions of the magnetic shear angle, and the asym¬ 
metric outflow speed and Sweet-Parker reconnection rate 
models. 

5. In simulations with a dipole tilt of +15°, we find most 
of the models map the magnetic separators when the IMF 
clock angle is 120°. However, none of the models chosen 
for this study faithfully reproduce the entire magnetic sep¬ 
arator for a northward IMF orientation of 30°. Again, only 
the maximum magnetic shear model traces the portion of 
the separator nightward of the magnetic nulls. The dayside 
portion of the magnetic separator in this simulation moves 
duskward by approximately 45° in longitude, as measured 
from the +®gsm axis. 

6. Accounting for flow shear parallel to the reconnect¬ 
ing components of the magnetosheath and magnetospheric 
magnetic fields improve the performance of the asymmetric 
outflow speed and asymmetric reconnection rate models un¬ 
der northward IMF orientations. Predictions for southward 
IMF orientations are mostly unaffected. 

In short, most of the models tested give results fairly close 
to the location of dayside reconnection for southward IMF 
orientations, but none are perfect. We conclude that the ex¬ 
isting models may be good enough for a ball park estimate, 
but it remains unclear what physically controls the location 
of magnetic reconnection at the dayside magnetopause. 

The present study employed a few underlying assump¬ 
tions. We employed a uniform, explicit resistivity in our 
global MHD simulations even though Earth’s magnetopause 
is collisionless; this choice ensures our simulations are well 
resolved, while reducing the likelihood of Flux Transfer 
Events (FTEs) [Russell and Elphic, 1978] in our simulations. 
However, recent advances have been made to trace magnetic 
separators in simulation when FTEs are present (A. Glocer 
et al., manuscript in preparation, 2014). 

Most of the reconnection location models discussed here 
are independent of dissipation mechanism, so the present 


results on the validity of each model should be independent 
of our choice to include an explicit resistivity. However, it 
is still unknown if the magnetic separators depend on dis¬ 
sipation mechanism. Future work will need to extend the 
present study to global Hall-MHD, hybrid, or kinetic simu¬ 
lations. 

The analysis presented here has largely focused on the 
global process of magnetic reconnection by identifying the 
magnetic separators across the dayside magnetopause. At 
present, there is debate as to whether magnetic reconnec¬ 
tion is a global or local phenomenon at the dayside mag¬ 
netopause. We determined a model’s prediction by maxi¬ 
mizing a parameter related to reconnection at the dayside 
magnetopause. Our underlying assumption is that a partic¬ 
ular model determines reconnection if that model faithfully 
reproduces the entire magnetic separator. Future work will 
need to assess this assumption by characterizing the loca¬ 
tions on the magnetic separator where magnetic reconnec¬ 
tion occurs. 

Finally, the present study uses a small set of solar wind 
and magnetospheric conditions. The present work has ig¬ 
nored the effect of including an IMF B x component. Previ¬ 
ous studies found that under southward IMF orientations, 
the reconnection site moves northward for B x > 0 and 
southward when B x < 0 [Peng et al., 2010; Hoilijoki et al., 
2014]. Additionally, we perform our analysis after the sim¬ 
ulations have achieved steady-state, which does not capture 
the magnetosphere’s response to dynamic solar wind con¬ 
ditions [Laitinen et al, 2006, 2007]. Understanding the re¬ 
sponse of Earth’s magnetosphere for a broader range of so¬ 
lar wind conditions is of the utmost importance for realistic 
space weather forecasting, and will be the subject of future 
work. 
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